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ABSTRACT 

As Space Shuttle tanks are loaded with N2O4, pressurization gas is 
displaced and at the same time N2O4 entering the tank evaporates and 
mixes with the pressurization gas remaining in the tank. Further addition 
of N2O4 to the tank requires venting of this mixture of N2O4/NO9 vapor 
and pressurization gas, which must be scrubbed prior to atmospheric 
release. 

A computer analysis was performed to estimate concentrations, flow 
rates, and total quantities vented during the planned fast fill/slow fill 
cycles of 125 gallon and 640 gallon Space Shuttle tanks. With helium as 
the pressurization gas, total vent quantities were estimated to be 4. 5 and 
24 pounds N 2 C>4, respectively, for the two tank sizes. 
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ESTIMATES OF SHUTTLE LOADING EMISSIONS 

by 

Thomas E. Bowman 
Mechanical Engineering Department 
Florida Institute of Technology 

1. STATEMENT OF THE PROBLEM 

The functional parameters pertinent to Space Shuttle hypergol loading 
operations are given in Table I. Essentially we are concerned with two 
types of tank for each propellant: 125 gallon tanks that are fast filled at 
35 psig and then slow filled at a pressure that rises to 125 psig and stays 
there, and 640 gallon tanks that are fast filled and slow filled at 35 psig. 


TABLE I 

Hypergol System Functional Parameters 


Orbital 

System 

Propellants On- 
Board (Gallons) 

Flowrate (GPM) 


Ullage Pressure 
(PSIG) 

MMH 

n 2 o 4 

Fast Fill 
0-90% 

Slow FiU 
90%-100% 

Fast 

Fill 

Slow 

FiU 

FRCS 

125 

125 

10 

2 

35 

125 

LH ARCS 

125 

125 

10 

2 

35 

125 

RH ARCS 

125 

125 

10 

2 

35 

125 

LH OM3 

640 

640 

60 

10 

35 

35 

RH OMS 

640 

640 

60 

10 

35 

35 

PBK 

3x640 

3x640 

60 

10 

35 

35 
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Assume that we start with a propellant tank full of helium at 35 psig, 
70°F, and begin to flow liquid propellant into the tank at a constant rate. 

As the first propellant enters, evaporation occurs, and liquid continues 
to evaporate as the tank fills until the gas over the propellant is saturated. 
At the same time, the vent is open and propellant vapor and pressurization 
gas are removed from the tank at a rate sufficient to maintain the desired 
ullage pressure. An estimate of the total quantity of propellant vapor dis- 
charged through the vent, as a function of time, is required. 


FIRST APPROXIMATION 

A very simple analysis may be performed by assuming that the mass of 
propellant vented is simply equal to the mass of saturated propellant vapor 
that will completely fill the tank minus the amount in the ullage at the end: 


m 


p. v. 


p sat ( V " V ull 
_____ 

P 



where V is the total tank volume, V u jj the final ullage volume, T the 
temperature of the gas mixture, and Rp and p ^ the gas constant and 

saturation pressure for the propellant in question. This approximation is 
based on a simplified model in which the tank is assumed to fill instantly 
with saturated vapor — without any venting being required — and then the 
saturated vapor is steadily pushed out the vent by the entering propellant. 

No further evaporation occurs because of the saturation condition. 

Saturation pressure of N 2 O 4 at 70°F is 14.71 psia, and the molecular 
weight (corrected for dissociation to NO 2 , assuming an equilibrium 
mixture of N 2 0 4 - NO 2 at 70°F, 35 psig) is 87.82, giving R p = 17.597 

ft-lbf/lbm - °R. We are concerned with two tank sizes: (V - V u u) = 125 
gal and 640 gal. The results are: 


125 Gallon N 2 0 4 Tank: 


m. 


p. v. 


= (14.71) (144) (125) (.1337) = 3.798 lbm 

(17.597) (529.67) 


640 Gallon N 2 0 4 Tank: 


I 
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-mp.v. = ^12_\ (3.798) = 19.45 lbm 

It should be noted that these calculations do not take account of the 
increased pressures during the last 10% of fill of the 125 gallon tanks. 

In the actual fill, the vent would close for a period of time to allow 
pressure to increase, and hence the quantity vented would be less than 
that calculated. After compression, the ullage gas still contains the same 
amount of propellant, but it is less than saturated because of increased 
temperature. No further evaporation occurs, however, because the liquid 
propellant is still at 70°F and its vapor pressure is now less than the 
partial pressure of the propellant in the compressed ullage. 

Beyond this consideration, however, it is not possible to ascertain 
whether this model is conservative or not. It has a conservative aspect, 
in that the gas being vented always contains the maximum possible con- 
centration of propellant. It also has an anti-conservative aspect, in that 
the gas flow rate out of the tank is always taken to be its lowest possible 
value — equal to the volume flow rate of liquid into the tank. Any time 
evaporation is occurring, the volume flow of gas out of the tank must ex- 
ceed the liquid volume flow in — a fact not taken into account by this 
always-saturated approximation. 

Because of the fact that evaporation of the propellant in the tank might 
cause the results found above to be substantially too low, a more careful 
analysis was undertaken. This analysis is best described as a "second 
approximation, " however, for reasons discussed later. It represents a 
preliminary estimate of the vent propellant quantities. 


3. SECOND APPROXIMATION 
3. 1 Formulation 

At any point in time, the volume occupied by gas in the 
tank is decreasing due to the addition of liquid to the tank, corrected 
by the amount of liquid that is evaporating: 



where 
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V(t>- 
Qin = 

•^evap 
liq = 


Tank ullage volume 

Volume flow rate of liquid 
into the tank 

Rate at which liquid is evaporating, 
mass per unit time 

Liquid density 


The volume flow of gas out the vent is equal to the rate at which the ullage 
volume is decreasing, plus a term to account for the generation of gas due to 
evaporation. We assume that the temperature and pressure of the gas mixture remain 
constant during evaporation, and that both components act as perfect gases. Consider 
a system composed of the gas molecules near the liquid-gas interface, with new 
molecules entering the system as liquid evaporates and no molecules leaving the 
system. This system will therefore be expanding. Its total volume can be treated 
as the sum of the partial volumes of each component gas because of the perfect gas 
assumption; the partial volume of the helium is constant (constant mass, temperature, 
and mixture pressure), while the partial volume of the propellant increases because 
of mass addition: 


V p = M p R p T/p 


ffi-T ^evap R p T/p 


where 


m F 

B P 

T 

P 


= partial volume of propellant vapor 
= mass of propellant vapor in the system 
= gas constant for the propellant vapor 
= mixture temperature 
= mixture pressure 


Since the partial volume of the helium is constant, the total volume of this system 
is increasing at a rate given by dV p /dt. This term is added to the rate at which 
the ullage volume is shrinking to give the volume flow of gas out the vent: 


Qout 


Rp T/p 


( 2 ) 


dV . 

_ dt~ m evap 

The portion of this volume flow that is propellant vapor varies with time 
because of the non-uniform distribution of propellant vapor in the ullage. Since 
the vent is the farthest point in the tank from the evaporating liquid surface, the 
propellant vapor concentration there should always be less than or equal to the average 
(uniform distribution) concentration, and hence a conservative result will be obtained 
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by assuming a uniform distribution of propellant vapor in the tank at all times. 

With this assumption, the density of propellant vapor at the vent is 

(3) 

RpT 

where the partial pressure of the propellant vapor, which the assumption says is 
uniform, is given by 

p = m R T/ t , (4) 

P p p V * ' 

nip being the total mass of propellant vapor contained in the ullage. The mass flow 

of propellant out the vent is the product of the total volume flow at the vent and the 

propellant vapor density at the vent: 


m r 


out 


= Qout 


(5) 


Finally, we have a continuity equation which provides a second relation 

between nu and m D : 
p out p 


dm, 

“dT 


p m evap “ m Pout 


( 6 ) 


Equations (1) through (6) constitute 6 equations - including first order 
differential equations - for the 6 unknown functions of time, V, Q out , _y^P> 

Pp, nip, and irip^, in terms of the parameters m e vapf Qin> y^liq » T »P» R p* 

Equations (2) through (6) are readily combined: 


dm P m evap ^out y* p 


l e vap “ | 

f- dV + 

k dt 

m 

9 

m evap R p 

t/ pJ 

( m p/V ) 


R T \ 

m 

dV 

evap ( 

'i- JL 

P ) 


K V 

P J 

V 

dt” 


(7) 


Equations (1) and (7) now constitute two first order differential equations for 
the two unknowns mp and V; for solution, two initial conditions will also be required. 

Before looking at the initial conditions, we should examine n^yap more closely. 
Evaporation is expected to be maximum when there is no propellant vapor. For the 
purpose of this analysis, then, we assumed m eva p to be given by an equation of the 
form 
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m evap A 


P , (T) 
sat 


- 1 i 


where 


P sa t ( T ) = saturation pressure for the propellant 
A = "Evaporation rate coefficient" defined by this equation 

Substituting this expression into equations (1) and (7), with p p given by (4), 

results in . 

/ \ 

V 


dV = 

Qin - _A_ 

dt 

' ^ liq 


/ P v \ 

dm = A 

/ sat -1 ] 


sat 


-1 


nipRpT 


( 8 ) 


1 - 


dt 


W 


m R T 
P P 

Tv - 


m ... 

+ _E. AX. 

V dt 


( 9 ) 


3.2 Initial Conditions 


Two initial conditions are required. One condition, V=V tank at 
t = 0, is quite straightforward. If we try to apply the condition mp =0, t=0, 
however, there is an obvious problem with equation (9). This problem can 
be resolved by considering that as t -*-0, mp ->-0 and p p ->0, and therefore 


p sat < T > 


r hevap 
so that with 


- A p sat (T) / (>H p T/ v tank 


Qhut = 0 at t = 0, we have: 


lim 
t 0 



A Psat Vtank 
m pBp T 


which can be readily integrated to give 


lim 
t -y0 


j m p) 

| 

^ 2A Psat Tank 

t 1/2 

(10) 


1 

R T 

'n 

■i 



which can be applied at some small finite time St, and hence is quite appropriate 
for the finite difference method of solution that will be applied to obtain the 
simultaneous solution of equations (8) and (9)* 

There are two limitations on the initial mp that must also be 
considered: 


67 



1. The initial m cannot exceed the total propellant mass 
that has entered the tank in the initial time step £t, 

m p 4 Q in Shq St (10') 

2. The initial nipCannot exceed the mass of propellant 

present when the tank is full of saturated vapor, 

P (T) V 

m P 4 ~f — — - — (10") 

P 

3.3 Slow Fill Analysis 

The slow fill presents no problems in the 640 - gallon tanks, 
where the fill and outflow rates simply decrease after 90% loading 
with no change in ullage gas properties. In the 125-gallon tanks, 
however, the vent closes for a period while the pressure increases. 

During the vent -closed period, we assume that the gas under- 
goes an isentropic compression, with no evaporation or condensation 
occurring. We further assume that it is a perfect gas with constant 
specific heats. These assumptions yield the equations 



for final volume and temperature, V 2 and T 2 , in terms of volume 
Vp temperature T-^, and pressure pj at the beginning of the 
compression, and final pressure p 2 . ~Y is the ratio of specific heats, 
which for a mixture of perfect gases can be expressed as 

<f 

j~ Pi C pi (13) 

p i 

= partial pressure of the i^ 1 component 

= the molal specific heat at constant pressure 
of the ith component 

= the molal specific heat at constant volume of 
the i*h component 


r= 


where 




Pi 


A 
C . 


VI 
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At standard temperature and low pressure, the specific heats 
have the following values: 

Helium C p = 5.00 Btu/lb mole - °R 

c v = 3.000 Btu/lb mole - °R 

N 2°4 1 c = (.2033) (87.82) cal/gmole -°K 

c“ v •= 'c - R = [(.2033) (87.82) - 1.986] cal/gmole -°K 

It should be noted that in the case of N 2 0 4 , the frozen equilibrium 
heat capacity was chosen as being the only one consistent with the perfect 
gas assumption. The equilibrium molecular weight at 70°F (87.82) was 
used in the equations above .to obtain the molal heat capacity from Fan 
and Mason's value of .2033 cal/gram - °C. 

Partial pressure of the N 2 0 4 vapor is given by equation (4) above. 
Helium partial pressure is found by subtracting the appropriate pro- 
pellant vapor partial pressure from the tank pressure. The resultant 
ratio of specific heats, equation (13), is: 

^N 2 0 4 = (.2033) (87.82) p N2 o 4 +5.00 (p -pj^ ) (13') 

tank 

ullage [(.2033) (87.82) - 1. 986]pj^o 4 +3, 000(p - Pj^O^.) 

If we use the saturation pressure at 70° F for the N 2 C >4 vapor partial 
pressure, we obtain the value 

^ N 2 0 4 = 1.293 

ullage 

and substituting this value into equations (11) and (12), with 
V 1 = 17.91 ft 3 - (.90) (125 gal) (.1337 ft 3 /gal) 

= 2.869 ft 3 

T x = 70° F = 529. 67°R 

Pj^ = 35 psig =49. 7 psia 

p 2 = 125 psig = 139. 7 psia 

gives the results: 

v - , . 1 

N 2°4 = 2 - , 4 ' J ' 7 1,393 290 ft 3 

ullage \ 139 - 7 / 

I 

1. Fan, Stephens. T. , and David M. Mason, "Properties of the System N 2 0 4 2N0 2 

2NO+ O9. " J. of Chem and Eng. Data, Volume 7, pp. 183-186 (April 1962). 
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3.4 


.293 
T 

2 

N 2 0 4 = 529.67 

ullage 

From these results we deduce two important facts: 

1. Since the final ullage volume, in the 125 gallon tanks is 

17.91 - (125) (.1337) = 1.20 ft 3 

maximum pressure will be reached before the tank is "full", and the 
vent will re-open to release some gas at 125 psig. 

2. The saturation pressure of an equilibrium mixture of and 

N0 2 calculated for 210°F is 14,920 mm Hg or 288.4 psia, and since 
the actual partial pressure of this constituent will increase in proportion 
to the tank pressure, the final partial pressure is well below saturation, 
and hence condensation at the relatively cool liquid and solid surfaces, 
although it will occur, can probably be neglected. 


139.7 X 1 - 293 


49.7 


= 669. 4°R = 210°F 


Numerical Solution 

Equations (8) and (9) were written as difference equations: 


A V = 

A m p = 


- >% 


P V 
sat 


Psat 


V 


m p R p T 


-)) 


.A t 


- 1 


m P R P T 


1 - m P R n T ) A t + 
pV 


(14) 


m p A>v 
V 


(15) 


and solved incrementally; that is, at each time step, the changes in V 
and mp were calculated using the values of V and nip calculated at the 
preceding time step. Initial values of V and Mp were found as de- 
scribed in Section 3.2. 

Values used were the following: 


P 

T 

At 


= 7156 lbf/ft 2 

= 529. 67°R 

= 1.0 sec 


/* liq = 


Psat 



90.35 lbm/ft 3 
2118.24 lbf/ft 2 
17.597 ft - Ibf/lbm - °R 
17.85 cal/mole - °K 
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Calculations were performed, in all cases, for five values 
of the "evaporation rate coefficient" defined on page 9: 


A = .01, .1, 1.0, 10.0, 100.0 lbm/sec 

(The effect of varying surface area on evaporation rates was not taken 
into account). The following sets of values were used for the two tank 


sizes: 

Initial 

v tank 

Final 

Propellant 
Volume (liq) 

• 

(^in)fast 

(Qin^ slow 

Small tank 

17.91 ft 3 

125 gal 

10 gal/min 

2 gal/min 

Large tank 

89.5 ft 3 

640 gal 

60 gal/min 

10 gal/min 


The procedure was as follows: 

1. For the initial time step, an average mass was calculated using equa- 
tion (10). The average was found in the usual way, 

(mp)ave = 


At 


At 


/ 


m (t) dt 

r 


_ __4_ f A Psat v tank A t \ 1//2 
3 ( R p T ; 

2. The limiting value of from eqn (10') was calculated and compared 
to the preceding value. If smaller, it replaced it. 

3. The limiting value of mp from eqn (10 M ) was calculated and compared 
to the preceding value. If smaller, it replaced it. 

4. Outflows for the first time step were calculated using the resultant 
value of mp. 

5. Volume and mass increments were calculated from equations (14) and 
(15) using these initial values of V and nip, and used to find new values 
of V and nip. 

6. This finite difference procedure was continued until either the ullage 
was saturated with propellant vapor, or the liquid volume in the tank 
reached 90% of its final value. Flow rates were calculated at each step 
using equations (1), (2) and (5), and a running total 


2- “Pout At 
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was kept. At each step, a check was made to ensure that the ullage 
had not become saturated and that a negative value was not being used 
for m eva p . 

7. As soon as the ullage became saturated with propellant vapor, the 
finite difference procedure was terminated. The remainder of the fast 
fill was analyzed using equations (1) through (G) with 

”Wap - 0 • 

In this case the equations are readily solved analytically, and numerical 
methods are not required. 

8. Following the completion of the fast fill, the slow fill was analyzed. 

For the 640 gallon tanks, this simply amounted to changing the value of 
Qj n and preceding as before. For the 125 gallon tanks, an isenti'opic 
portion was calculated using equations (11), (12), (13a) and (13b), after 
which the last stage of venting was analyzed as described in steps 6 and 7, 
with p . = 20116 lbm/ft^ 

T = Tg (just calculated) 

Q in = 2 gal/m in 

and starting with V = V 2 , the value found for the isentropic compression. 

The results for total propellant vent masses are summarized in 
Figures 1 and 2. More detailed results, including a listing of the computer 
program, are presented in the Appendix. 

4. Discussion of Results 

No attempt was made to select realistic values of A, which will depend 
on the liquid propellant surface area and hence will vary with time, especially in 
the case of a spherical tank. The results varied with A throughout the range of 
values of A that were investigated. It is likely that realistic values of A will 
fall within this range; however, a more thorough study would be required to de- 
fine the vent masses with a high degree of certainty. 

Two other shortcomings of the analysis should be pointed out: 

1. Some discrepancies exist in the way V, nip, and m ga t are handled 

in the first time step. . 

2. The slow fill analysis for the 640 gallon tanks does not allow for 
the possibility that the ullage might not be saturated. 
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It should be emphasized that this analysis was intended only as a "second approxi- 
mation, " and was a fast-reaction exercise of very limited duration. We did not have 
time or funding, at the time this work was done, to refine the analysis beyond the 
point where consistent and reasonably reliable results were being obtained. Our 
approach and results. Including emphasis on the deficiencies of the analysis, are pre- 
sented here since we are not aware of the existence of anything better, and for the 
benefit of anyone who might wish to develop this computer code further. 

Figures 1 and 2 are cross-plots of the total amounts vented for each tank versus A. 
Figure 1 in particular indicates that the propellant mass vented might well grow without 
bound as the evaporation rate increases. At these high values of propellant vapor mass 
however, it would be necessary to take account of cooling due to evaporation; this cooling 
would reduce the vapor pressure and hence the amount evaporated. In addition, we should 
remember that these high values occur as a result of the very rapid evaporation at early 
values of time, accompanied by an assumed instantaneous diffusion throughout the tank, 
with the result that large amounts of propellant are lost out the vent before the gas is 
saturated. A final consideration: examination of the computer output shows that the curves 
would stop rising at only slightly higher values of A than those calculated because of the 
limitations represented by equations (10') and (10 M ). 

With these considerations in mind, it was decided to select the approximate inflection 
points in the two curves to estimate the total vent quantities; 

125 gallon N 2 0 4 tank: m p> v . = 4. 5 lbm 

640 gallon IS^O^ tank: m p y = 24.0 lbm 

The 640 gallon figure corresponds to a higher value of A than the 125 gallon figure, 
which is appropriate since A should be proportional to surface area. 

These values are 18.5% and 23.4% greater, respectively, than those obtained with 
the first approximation. They are also 5 - 10% smaller than the maximum values 
calculated. 

The complete flowrate history (hypergol vapor vent flow rate versus time) calculated 
by this program is presented in Figure 3. 
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TOTAL PROPELLANT VENTED (LBM) 
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Figure 3: Baseline Flowrate History, - NO^ vapors vented during pad servicing operations - does 

not include helium flow 








APPENDIX 

VENT FLOW RATE PREDICTIONS FOR LC-39 
(NORMAL OPERATING CONDITIONS ) 


This Appendix includes a listing of the computer program used 
for these calculations, and results for selected values of the 
"evaporation rate coefficient, " A, defined in Section 3. 1. 
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DIMENSION A<5).V0(4) . RHG < 2 ) . PS AT I 2 > . R ( 2 ) . T I TLE ( 1 0 > . QDOT Jt 4 ) . 
1 VTOT<4) , CO( 4 1 
I N=5 
I OUT-6 

READ ( IN, l ) IA( I).I=1 ,5>.( V0( J),J=1.2> 

READ < I N , 3 ) I VTOT ( J ) . J= 1 .2 t 
READ ( I N • 3 ) ICDOTJt J) ,'Js| ,4) 

PFAO ( IN , 3 I ICP(K).K= I. 2) 

3 FORMAT ( 8F 1 0 • 3 > 

T — 529, 67 
PM I X= 7 156. 

PMIX2=201 16 
DELT=1 . 

RHO ( I >=54. S 
RH0(2)=90.35 
PSAT< 1 1=1 1 1 .4S6 
PSAT( 21=21 1 e.24 
R(l)=33.536 
R(2)=t7.S97 
DO 100 K=1.2 

REAO I IN. 21 (TITLE(L) .L=1.9) 

WRITE (IOUT, 10) (TITLE(J).J=1.9> 

DC 100 J=1 .2 

WRITE { l OUT ,11) VO ( J) 

DC 100 1=1 ,S 

WRITE C IOUT, 12) A( I 1 

WRITE 1 1 OUT ,14) 

T I ME = DEL T 

KCUN T= 1 . , • 

A A = A ( I 1 

ULLAGE = VTOT { J )!— . 9WV0(J)*.1337 
OOOT=ODOTJ ( J) 1.337/60 . 

V= VTOT I J 1 

EMM=CDCT*DELT*RHO( K) 

EMM2=S0RT(a./9.*AA*PSAT{K)*V*DELT/R<K)/T) 

IF (EMM2.LT.ERMI EMR=EMR2 
EMM?=PSAT (K)*V/R(K)/T 
IF (EMM2.LT. EMM) E«R=EMM2 

FLOW=ODOT+EMM*( R(K) 4T/PM IX-l ,/RHCIK) l/DELT 
FLAW=FL0W/ . 1 337 

EFFLUX=FLCW*EMM/V . 

SCRU9=EFFLUX*DELT 
44 PARA=PSAT (K > * V/EMM/R ( K )/T 
IF (PARA. LE. 1.1 GO TO 66 
E VAP= A A * ( P AR A- 1 . ) 

TEST=0D0T*RHG(K) 

IF (EVAP.GT. TEST) EVAP=TEST 
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OELV=<EVAP/PHC(K)-CCnT)*OELT 

DELM=EVAP*< 1 . -PSAT (Kl/PMIX/PARA J*CELT+EMM/V*OELV 

22 TIME = T1ME+QF.LT 

t WM=E MM+DELM 

V=V+DELV 

PP=EMU*P(kI»T/V 

FLOVi=QDf)T + EVAPA{P{K)*T/P«lX-l ./BHC(K> > 

FLAW=FLOW/. 1337 

EFFLUX=FLQW*EMM/V 

SCRUB — SCPUH* EFFLUX *DE LT 
KCUNT=KCUNT+ l 

IF <KOUNT.LT. 101 GO TO 33 

VOUT=V/.1337 
POUT — PP/144. 

WRITE ( I OUT , 1 3 > T I ME . E MM. V CUT . POUT. FLAW, EFFLUX. SCRUB 
KCuNT _ 0 : ; 

33 IF (V.GT. ULLAGE > GO TO 44 
VOUT— V/« 1337 
P0UT=PP/144. 

GO TO 88 

66 EMM=PSAT(KI*V/R(K)/T 

_________ ■ — - — — r • 

VCUT=V/.1337 

P0UT=PSAT(K)/1A4. 

" write T I out . 13) T I me, eMM, v(3ut ;PCuT;PLWrEFFLUX^SCRiJ6 ■ 

CHEKl=PSAT(K> *ULLAGE/fi<K)/T 
TIOE=( V-ULL AGEI/OOOT 

T ime^timEatioe ' ' “ ' ~ ' 

VOUT=ULL AGE/ . 1337 
P0UT=PSAT(K)/I44. 

a W = 0'00T/V 1 3 5 7 * 

EFFLUX=aOCT*PSATIK)/R<K»/T 
SCRUP=SCPun*EFFLUX*T IOE 

C HE K2=E M M— EFFLUX *T ICE I ~ 

IF I ABSICHEK I-CHEK2 > .GT. .01 4CHEKI ) WRITE (I0UT.15) 

£MM=CHEK1 

88 wp i te - < Tout , i3> t ime.emm.vcut.pcut, flaw. efflux, scrub 
OOOT=QDOT J ( J* 2 ) * • 1337/60. 

ULOLD=ULLAGE 

ULLAGE* VTOt ( J ) -VO { J I 4 . I 33~7 "" " 

IF ( J.EQ. It.GC TC « 

EMK=PSAT< K) AULLAGE/PIK J/T 
T IDE* VO I J > *.0 1337/QCOT 
TIW6=TIFE+TI0E 
VCUT=ULL AGE/ • 1337 
FLAW=ODOT/. 1337 
EFFLUX=COOT*PS AT { K) /P <K)/T 
SCRUB* SC RU8 A EPF LUX* TICE 

WRITE ( I OUT. 13) TIME, EMM, VOUT.POUT, FLAW, EFFLUX, SCRUB 
GO TO ICO 
99 WRITE IinUT.lE) 

GAMMA=(CP*CP(K)+ (PM IX-PP )*5«)/<PP*(CP<K»-l.<SS6J + (PMlX-PpJ*3.» 
PCWI =1 ./GAMMA 

V I SEN=ULOLO* ( PMI X/PM 1X2 ) **PCW 1 
IF ( VISEJN.LT. ULLAGE) GC TO TOO 
T I ME = T I ME* I ULCLO- V ISEM/GOCT 
VOUT=VISEN/. 1337 


PP-PP* PMI X 2 / P M I X 






POUT— PP/1 44 . 

FLAW=0. 

, EFFLUX=Q. ■ ’ . 

J ■ WRITE (IOUT. I 3 ) TI M E. EMM, V OUT, POUT, FLAW. EFFLUX, SCRUR 
WRITE ( 1 CUT ,17) 

• POW2= (GAMMA- I . ) /GAMMA 

T2=T*(PMI X2/PM IX)**POw£ " 

TIOE=( VI SEN-ULLAGE I /OCOT 
T IME=TIME + TIOE 
VOUT=ULL AGE/ . 1337 
FL AW=OCOT/ . 1 337 

efflux=coot*pp/r < k ) /t 2 

SCRUS=SCRUR*CFFLUX*TIOE 

E MMsEMM— EFFL UX4T I CE 

WRITE (IOUT. 13) TIME. EMM, VOUT.PCUT, FLAW, EFFLUX. SCRUB 
100 CONTINUE 
STOP 

1 FORMAT (SF10.3/4F10.3) ' ■ 

2 FORMAT (9«4) 

10 FORMAT ( 1 H I ,9 A4 ) 

11 FORMAT {///• FINAL PROPELLANT VOLUME =',F6.0,* GALLONS*/) 

12 FORMAT (/■ EVAPORATION RATE COEFFICIENT =*,F7.2, *LBM/SEC*/> 

13 FORMAT (FB. 0 ,F I 1 . 3, FI 1 .2.F l.C .3. FI 1 .3. FI 2.4 ,F12. 3> 

1A FORMAT ( //3X , • T I ME PROPELLANT ULLAGE PARTIAL TOTAL PR 

1R0PELLANT* ,6*,< TOTAL* /3 X , • (SEC) VAPOR MASS VOLUME PRESSURE 
2 OUTFLOW OUTFLOW PROPELLANTS •/ 1 AX .*( LRM ) (GAL) (PSIA 

3A> (GAL/SEC) (LEM/SEC) VE NT EO ( LBM ) • / ) 

15 FORMAT ( • OOPS* I 

16 FORMAT (• VENT CLOSES') 

17 FORMAT (• VENT OPENS' ) ' ■ : “ “ ~ ’ ' 

END 


80 


N 2 0 4 RESULTS - 125 GALLON TANK 


EVAPORATION RATE COEFFICIENT ' = 1.00LBM/SEC 


TIME 

(SEC) 

propellant 
VAPOP MASS 
(LBV) 

ULL AGE 
VOLUME 
(GAL) 

PART ! AL 
PRESSURE 
(PSIA) 

TOTAL 
CUTFLCW 
( GAL/SEC) 

PRCPELLANT 

OUTFLOW 

(LBM/SEC) 

■ fofAL 

PROPELLANTS 
VENT. ED I LBM ) 

10. 

3. 862 

132.66 

14.095 

0.691 

0.0201 

0.876 

20. 

3.958 

131.01 

14.626 

0.2 34 

0. 0071 

0.984 

30. 

3.927 

129.35 

14.698 

0.1 76 


1 .043 

AO. 

3.879 

127.68 

1 4. 708 

0.1 68 

0.0051 

1.094 

50. 

3.829 

126.01 

14.710 

0.167 

0.0051 

1.145 

se. 

3.788 

12A.68 

14.710 

0.167 

0.0051 

1.186 

677. 

0.652 

21.46 

14.710 

0.167 

0.0051 

4.322 

VENT CLOSES 

1032. 0.652 

9.65 

41.361 

0.000 

0 .0000 

4.322 

VENT OPEN 
1052. 

3 

0.605 

8.96 

41. 351 

0.033 

0.0023 

4.369 
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N 2 0 4 RESULTS 640 GALLON TANK 


EVAPORATION PATE COEFFICIENT = 10.00LBV/SEC 



TIME 

PROPELLANT 

ULL AGE 

PART I AL 

TOTAL 

PRCP6LL ANT 

TOTAL 

( SEC ) 

VAPOR HASS 

VOLUME 

PRESSURE 

OUTFLOW 

OUTFLOW 

PROPELLANTS 


(LBN) 

(GAL 1 

(PSIA) 

{ CAL/SEC > 

(LBM/SEC) 

VENTED ( LBM) 

10. 

20.029 

661.26 

14.663 

1.477 

0.0447 

5.1 12 

20. 

19.786 

651.27 

14.709 

1 . 0C6 

0.0306 

5.441 

30. 

19.485 

641.27 

14.710 

1.000 

0.0304 

5.746 

31 . 

19.455 

640.27 

14.710 

i .ooo 

0.0304 

5.7 76 

578. 

2.838 

93.4 1 

1 4. 71 0 

1 .000 

-0.0304 

22.392 

96 2. 

0.894 

29.41 

14.710 

0.167 

0.0651 

2T7337 


I 


Slow fill commences at 578 sec, ends at 962 sec. 
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1.000 2.000 
PROPELLANT OUTFLOW (LBM/SEC) 


lbm 

sec 



640 GAL. TANK, N^, EVAPORATION RATE COEFFICIENT - 10.0 





30.0 


4 25.0 


4 20.0 


H 15.0 


1 10.0 


H 5.0 


TIME(SEC) 


TOTAL PROPELLANTS 













